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Abstract: We propose a novel method to improve the spatial resolution of Brillouin optical 
time-domain reflectometry (BOTDR), referred to as synthetic BOTDR (S-BOTDR), and 
experimentally verify the resolution improvements. Due to the uncertainty relation between 
position and frequency, the spatial resolution of a conventional BOTDR system has been 
limited to about one meter. In S-BOTDR, a synthetic spectrum is obtained by combining 
four Brillouin spectrums measured with different composite pump lights and different 
composite low-pass filters. We mathematically show that the resolution limit, in principle, 
for conventional BOTDR can be surpassed by S-BOTDR and experimentally prove that 
S-BOTDR attained a 10-cm spatial resolution. To the best of our knowledge, this has never 
been achieved or reported. 

Keywords: Brillouin scattering; distributed sensor; optical fiber sensor; optical time-domain 
reflectometry (OTDR); strain measurement; synthetic spectrum; temperature measurement 



1. Introduction 

The spectrum of Brillouin scattering in optical fiber shifts in proportion to changes in the strain 
and temperature of the fiber. Brillouin optical time-domain analysis (BOTDA) [1,2] and Brillouin 
optical time-domain reflectometry (BOTDR) [3,4] are distributed measurement techniques utilizing 
this property by injecting a pulsed pump light and observing the scattered light. The amount of 
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strain/temperature change is estimated by measuring the spectral shift, while its position is determined 
by the round-trip time of the light. Since BOTDR uses only one end of a fiber, it is suitable for long 
distance measurement, whereas BOTDA uses pump and probe lights that are injected from both ends of 
a fiber. 

To improve the spatial resolution of both techniques, it is necessary to narrow the pulse; however, 
the linewidth of the observed spectrum then becomes wider and makes measurement of the spectral 
shift difficult. For this reason, the spatial resolution has been limited to about one meter in both 
techniques [5,6]. 

For BOTDA, Bao et al. [7] found experimentally a phenomenon in which the spectral linewidth 
becomes shorter when a very short pulse of about 1 ns is used. This phenomenon arises only when 
there is light leakage. Inspired by this discovery, various resolution improvement methods have been 
proposed [8-14]. The pump light of BOTDA plays two roles: phonon excitation and scattering by the 
phonons. The width of the pump light must be longer than the phonon lifetime (about 9 ns) for phonon 
excitation, whereas it must be much shorter than 10 ns, which corresponds to a one-meter resolution, 
for high spatial resolution. To satisfy these incompatible requirements, an idea to construct a pump light 
using a long and a short element was born. The long element, which is a long pulse or a CWwave, takes 
the role of phonon excitation, and the short element takes the role of being scattered by the phonons. 
A part of the spectrum is by the combination of these elements. By emphasizing this desired part, 
high-resolution measurement could be attained. 

For the construction of pump light from two elements, the authors of [8-11] used amplitude 
modulation and [12,14] used phase modulation. In any construction method, the measured spectrum 
includes both desired and undesired parts. Parameters must be optimized to emphasize only the desired 
part. Instead of that, the authors of [12-14] proposed methods to suppress or cancel the undesired 
part of the spectrum by combining two different measurements. By these methods, high-resolution 
measurement of centimeter-order for BOTDA has been attained. 

On the other hand, the pump light of BOTDR plays a single role, since the BOTDR utilizes 
spontaneous scattering, so that it is difficult to attain high resolution by the same idea as for BOTDA. For 
resolution improvement of BOTDR, Koyamada et al. [15] proposed a double pulse method that yields 
an oscillatory Brillouin spectrum and verified 20-cm spatial resolution by experiment. 

In this paper, we propose a novel method, referred to as synthetic BOTDR (S-BOTDR), to improve the 
spatial resolution of BOTDR. In this method, a synthetic Brillouin spectrum is constructed by combining 
several spectrums obtained by BOTDR measurements with different pump lights and low-pass filters. 
The pump lights and low-pass filters are composed of short and long elements with phase differences. 
We mathematically show that the resolution limit in principle for conventional BOTDR could be 
overcome by S-BOTDR and experimentally prove that 10-cm spatial resolution, which is much smaller 
than the conventional BOTDR limit, is attained by S-BOTDR. 

The remainder of this paper is organized as follows: Section 2 presents the mathematical formulation 
of BOTDR and its performance limit in principle. In Section 3, the proposed method, called S-BOTDR, 
is described. Sections 4 and 5 are devoted to the evaluation of the proposed method through simulations 
and experiments, respectively. Finally, we conclude our paper in Section 6. 
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2. Mathematical Model of BOTDR 

2.1. Sensing Mechanism of BOTDR 

A BOTDR system is shown in Figure 1. The light from the light source is divided into a pump 
pulse and a reference continuous wave. In an optical fiber, the acoustic waves are always excited by 
thermal fluctuation and the pump pulse is scattered by the acoustic waves in all parts of the fiber. Since 
the acoustic waves move with the sound velocity of the medium, the backscattered light has a Doppler 
shift, which is called a Brillouin frequency shift (BFS). As the sound velocity changes in proportion to 
the change in strain or temperature in each position of the fiber, the BFS also changes in proportion to 
them. In BOTDR, the backscattered light is heterodyned with the reference wave and the spectrum of 
the interfering light; i.e., the Brillouin spectrum is obtained. The BFS is estimated from the Brillouin 
spectrum as its center frequency. Since the round-trip time of the light corresponds to the position in the 
fiber, the BFS at each position in the fiber is obtained. Thus, distributed sensing of strain and temperature 
is possible using the Brillouin scattered light. 

Figure 1. Block diagram of Brillouin optical time-domain reflectometry (BOTDR). 
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2.2. Basic Equations of BOTDR 



Brillouin scattering in optical fiber is described by the following equations [16-18]: 



^ ^ + ^ + a 1 E 

Vo- dt dz 2 



Id d a . ^. 

^Vgdt dz 2 

?£+(^- + 2mvB{z))p = iAE p E* s +R(z,t) 



(1) 

(2) 
(3) 



where E p , E s and p are complex amplitudes of a pump light, a backscattered light and an acoustic 
wave, respectively, k and A are coupling coefficients, v g is the light speed in the fiber, a is the fiber loss 
factor and superscript * denotes a complex conjugate. T B is the linewidth of the acoustic wave spectrum 
and ^b(^) denotes the BFS at position z. Here, the spectral linewidth is defined by the full width at 
half maximum. 
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Random force R(z, t) is a circular symmetric complex white noise that is white both in space and time; 
i.e., it is characterized by E [R{z, t)R*(z', t')] = QS(z— z f )8(t— t') 9 where E[-] stands for expectation and 
Q is a constant. The two terms in the right-hand side (RHS) of Equation (3) correspond to stimulated 
and spontaneous Brillouin scattering, respectively. Under usual BOTDR conditions, as analyzed in a 
previous paper [19], the largest part of the measured spectrum arises from spontaneous scattering, and 
the stimulated scattering term can be neglected. Moreover, then, since it becomes \E S ] <C |-E7 P |, the RHS 
of Equation (1) can be neglected. Therefore, BOTDR is described by the following equations: 

1S +^ + ?U = 0 (4) 



v g dt dz 2 

^ + (T + 2mu B (z))p = R(z,t) (6) 



The boundary conditions of E p (z, t) and E s (z, t) are given by: 



E P (M = x /^f(t) (7) 
E s [ z, -) = 0 (8) 



where P p and f(t) denote the power and shape function of the pump pulse injected into an optical fiber, 
respectively, A e ^ is the effective core area of a fiber and T = T B /2 is set. 

2.3. Analytical Solution to the BOTDR Equations 

The solution to last section's BOTDR equations can be represented analytically. For simplicity, 
assuming that the fiber loss is small, we set a = 0. 

First, the solution to Equation (1) under boundary Condition Equation (7) is represented as: 



Next, the stationary solution to Equation (6) is represented as: 

p(z,t)= I e- {r+2 ™ Biz)){t - s) R(z,s)ds (10) 

J — OO 

whose autocorrelation function is given by: 

E \p(z,t)p*(z?,lf)] = Q;6(z - z')e~ 2 ™ B(*)(t-t') e -r|t-t'| 
Then, substituting (9) and (10) into (5) and solving it under Condition (8), we obtain: 

E s (z, t) = ik, r fit- 2 -^-\ p* ( z ',t- ^— ^ dz' (12) 



where ki = ^/P p /A eff /^ and L f is the length of the fiber. 
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The backscattered light returned to the input end of an optical fiber in BOTDR is represented as: 

X(t) d ^ E s (0,t)=iK 1 £ f f(t-^)C(z,t)dz (13) 
where we set ((z, t) = p(z, t — z/v g ), which has the same statistical property as p(z, t); i.e., 

E [C(*, t)C(z', t')\ = Q-5(z - z ') e -^B(z)(t-t') e -r\t-t'\ (14) 

holds. We note that X(t) becomes a circular complex Gaussian (ccG) process with mean zero. The 
component of this signal with frequency v is obtained by: 

Y(t,v) = c Y h(t)*[X(t)e- 2niiyt } 

= m 2 J h(t - r)e- 2 ™ T * f (r - ^ r)dzdr (15) 

where c Y is a constant, h(t) is an impulse response of a low-pass filter and * denotes a convolution, and 
we set K, 2 = cyKi. y(t, v) also becomes a ccG process with a mean of zero as X(t). 

2.4. Brillouin Spectrum and Point Spread Function 

The Brillouin spectrum obtained by one measurement of BOTDR is represented as: 

V(t,is) = \Y(t,is)\ 2 (16) 

Since BOTDR's signal source is thermally excited acoustic waves, its spectrum fluctuates even without 
observation noise. Although the spectrum is averaged by many replicate measurements, we must 
suppress the fluctuation for one measurement. 

Since Y(t, v) is a ccG process, V(t, v) becomes a random variable with an exponential distribution 
for each t and v, so its variance equals the square of its expectation: 

E [V(t, v) - EV(t, is)] 2 = (EV(t, v)f (17) 

To obtain a smooth spectrum, measurements are performed many times and the spectrums are 
accumulated or averaged. Therefore, the Brillouin spectrum obtained by measurements is considered 
to be expectation EV(t, v). 

We can calculate the expectation of V(t, v) by using the statistical property of Equation (1 1) as: 

^(t,i/)=7feL(t,i/) t * (18) 

where 7 R = v^ 2 Qj2Y 2 is a constant and ** denotes a two-dimensional convolution with respect to t 
and v. Functions Lit, v) and ^(t, v) are defined by: 

r 2 

LM = ri + (2n(v-Mv g t/2)y (19) 

*(*,^) = \F T [f{r)h{t-T)]\ 2 (20) 
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respectively, where T T denotes the Fourier transform with respect to r. Equation (18) is derived in 
Appendix A. 

We note that L{t,v) is a time- varying Lorentzian spectrum and is determined only by the 
characteristics of an optical fiber, whereas \^(£, v) depends on the sensing mechanism of BOTDR. Since 
Equation (18) implies that function, v), obscures the details of L(i, v), we refer to ^(i, v) as a point 
spread function (PSF). 

2.5. Performance Limit of BOTDR 

Though the ideal Brillouin spectrum is the Lorentzian spectrum, the observed Brillouin spectrum is 
spread by the PSF. Therefore, it is desirable for the PSF to have a shape close to a two-dimensional 
5-function; that is, narrow in both time and frequency directions. However, by the following uncertainty 
relation (see Appendix B): 

2A T • 2A F > - (21) 

TT 

the product of time and frequency widths of the PSF cannot be reduced below a certain value. Since 
the time width corresponds to the spatial resolution, there is a limit, in principle, to improving both 
spatial and frequency resolution simultaneously through one measurement with whatever pump pulses 
and low-pass filters are used. 

3. BOTDR by the Synthetic Approach 

To overcome the resolution limitation of conventional BOTDR described in the previous section, a 
synthetic spectrum approach, which we refer to as synthetic BOTDR (S-BOTDR), is proposed. 

3.1. Elements of the Point Spread Function 

We consider a short pulse, and a long pulse, f2(t), as elements of composite pump light: 

h{t) = I [t0 , t0+Dl] (t) (22) 
f 2 (t) = rI l0>Da] {t) (23) 

where is the definition function of interval [a, 6], Di and D 2 are the pulse widths of two pulses and 
r is the amplitude ratio of the two pulses. The start time of the short pulse is denoted by t 0 . 

Similar to the PPP-BOTDA method [8,9], Di is set as less than the spatial resolution to be attained 
and D 2 is taken as sufficiently longer than the acoustic wave lifetime, 2/T B , say about 9 ns, to improve 
frequency resolution. A low-pass filter is also composed of two elements, each of which is a matched 
filter of short or long pulses: 



hi(t) = h(-t) 
h 2 (t) = f 2 (-t) 



(24) 
(25) 
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Figure 2. Pump lights and matched filters: (a) Element 1 of the pump light; (b) Element 2 of 
the pump light; (c) Element 1 of the low-pass filter; and (d) Element 2 of the low-pass filter. 
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By using the above elements, a pump light and a low-pass filter are composed as: 

fe(t) = fi(t) + e id f 2 (t) 
K(t) = h^t) + e i4 >h 2 (t) 



(26) 
(27) 



where 9 and 0 are the phase differences of the two elements. The PSF corresponding to the composite 
pump light and low-pass filter becomes: 

*eAt,v) = |-M/*(r)M* - ^)}| 2 

= \F n \ 2 + \F 12 \ 2 + \F 21 \ 2 + |F 22 | 2 + [e~ i6 '(F 11 F 2 * 1 + F 12 F* 2 ) 
+ e-^(F n F* 2 + F 21 F 2 * 2 ) + e-^-^FuF^ + e' 1 ^ F n F* 2 ] 

where ] denotes the real part, and the RHS is composed of functions: 

F w = F kl {t,v)=F T [f k {r)h l {t-T)} 
= FAfkWfiir-t)], k,l = l,2 



(28) 



(29) 



The spreadness of function F k i has the following features, due to the difference in pulse widths of the 
two pulses (see also Figure 3): 

Fn: small in time and large in frequency, 
Fu and F 2 i. large in both time and frequency, 
F 22 : large in time and small in frequency. 

Therefore, among the terms of the RHS of Equation (28), only the term, F n F 2 * 2 , has the ideal property 
that it is localized both in time and frequency, and the other terms are undesired elements for the BOTDR 
measurement. 
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Figure 3. Illustration of F k i, k, I = 1,2 spreadness. |F fc ;|'s take large values in the black 
regions: (a) F n ; (b) F 12 ; (c) F 21 ; and (d) F 22 . 
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3.2. PSF by the Synthetic Approach 

To extract only element 5R(FnF 2 * 2 ) from the PSF, we prepare p phase pairs ((0j , j = 1, 2, • • • ,p) 
and p pairs of pump lights and low-pass filters: 



f (j) (t) = f!(t) + e»>Mt) 
h {j \t) = h 1 (t) + e i ^h 2 (t) 
J = 1,2,-.- ,p 



(30) 
(31) 



The Brillouin spectrum and the PSF obtained by using f^it) and h^\t) are denoted by V^(i, ^) and 
^\t, v), respectively. The synthetic spectrum is defined as: 



Vs(t,v) = Y i c j V<*\t,i,) 
i=i 

where Cj,j = 1, 2, • • • , p are weighting coefficients. The corresponding synthetic PSF becomes: 



(32) 



i=i 



(33) 



The j-th PSF is represented as: 



i&U\t,i>) = \7 T [f^{r)h {j \t-r 



(34) 
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which is expanded as: 



|^ll| 2 + |^12| 2 + |^2l| 2 + |^22| 2 

+ 23? [e-^(F u F^ + F l2 F* 22 ) + e-^(F n F* 2 + F 2l F* 22 ) 



If the following relation holds by combining the PSFs, an ideal PSF is obtained. 

^ s (t,u) = 2p^[F u F; 2 } 
This problem is equivalent to finding 9j, <fij, Cj, j = 1, 2, • • • ,p that satisfy: 

1 ••• 1 \ 



(35) 



(36) 



/ 1 

e -i6 



-upl 



-tt>2 



-Vp2 



e -i(dl-4>l) e -i(d 2 -<j>2) 



\ 



-i(0i+0i) 



-i(6>2+0 2 ) 



1 

e -i(0 p -cf)p) 
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( cA 
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/o\ 

0 
0 
0 

\p / 



(37) 



Though the general solution has not been obtained, the following solution with p = 4 is at least a solution 
to Equation (37): 



/ ft \ 

03 
V 04 / 
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(38) 



Hereafter, we use this solution. 

The PSF by the synthetic approach is represented as: 

*s(t,v) = 8U[F n F, 



22\ 



(39) 



and is concentrated around the origin in both the time and frequency directions, as shown in the previous 
section. However, in addition to this, its two-dimensional integral in the time-frequency space must be 
a certain positive value to approximate a two-dimensional <5-function. The two-dimensional integral of 
the synthetic PSF of Equation (39) becomes: 



/OO POO / POO 

/ * s (t,v)dvdt = 8[ fi{t)h{t)dt 
-oo J — oo \J —oo 



(40) 



If two pulses are separated in time, this integral becomes zero, and if they overlap partially, the integral 
is proportional to the width of the overlapping parts. To maximize the integral, the short pulse must be 
included in the long pulse, and the integral becomes 8r 2 D 2 . 

For two pulses of Equations (22) and (23), the synthetic PSF becomes: 



Vs(t,v) = 8r 2 cos7r(L> 2 - D r - 2t 0 )is 



sinvT^i — \t\)u sin tt(D 2 — \t\)v 



7T1S 



7TU 



I[-D lt D!](t) 



(41) 
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If Di is sufficiently small and D 2 is sufficiently large, the PSF is approximated as: 

tfsM ~ %r 2 D\5{t)5{y) (42) 
and the synthetic spectrum is approximated as: 

EV s (t,v)n8r*D 2 1 >Y R L(t,v) (43) 

The PSFs for a conventional BOTDR and S-BOTDR are plotted in Figure 4. While the PSF of the 
conventional BOTDR spreads in either the time or frequency direction, that of S-BOTDR is concentrated 
at the origin and approximates a 5-function. To clarify the mechanism of generating such a 5-like 
function, the composition of the synthetic PSF is shown in Figure 5. Although each of the PSF, 
= 1,2,3,4, concentrates in the time direction, it is spread in the frequency direction. This 
spreading is partially suppressed with two-by-two combination, \E^ 2 ) — and *( 4 ) — \E^ 3 \ and is 
completely suppressed with the use of all PFS elements. 

Figure 4. Point spread functions: (a) conventional BOTDR (D = 10 ns); (b) conventional 
BOTDR (D = 1 ns); (c) synthetic (S)-BOTDR (D 1 = 1 ns; and D 2 = 50 ns). 




3.3. S-BOTDR System 

As shown in the analysis of the previous section, the synthetic Brillouin spectrum has the ideal 
property of being close to a Lorentzian spectrum. An S-BOTDR system for obtaining the synthetic 
Brillouin spectrum is constructed as shown in Figure 6. The composition of the synthetic Brillouin 
spectrum is followed by Equation (32) and is illustrated in Figure 7. Although the spectrum of each 
measurement exhibits high spatial resolution, it is considerably different from the Lorentzian spectrum. 
Only the combination of all four spectrums produces the Lorentzian spectrum. 

The ideal property described in the previous section is for the expectation of the synthetic Brillouin 
spectrum. Achieving this property requires averaging or summing of the spectrums obtained by a number 
of measurements. 
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Figure 7. Composition of the synthetic spectrum of S-BOTDR. 




3.4. Parameter Selection for S-BOTDR 

The pump lights of S-BOTDR are composed of short and long pulses and include parameters t 0 and 
r: t 0 is the difference of pulse start times and r is the amplitude ratio of the two pulses. These parameters 
affect the performance of S-BOTDR. 

Selection of t 0 If the short pulse is included in the long pulse, i.e., if 0 < t 0 < D 2 — Du the integral 
of Equation (40) takes the maximum value. In addition, the value of t 0 affects the concentration of the 
PSF near the origin. The degree of concentration is estimated by the integral: 

/V POO 
/ u')dtdu'/(8r 2 Dl) (44) 

-v J —oo 

If C{y\ t 0 ) approaches one for small v, the PSF is considered to be concentrated near the origin of the 
frequency axis. 

The values of C{y\ to) are plotted for different t 0 's in Figure 8. Here, for £ 0 = (D 2 — Di)/2, the short 
pulse is located at the center of the long pulse; for t 0 = 0, the short pulse is at the edge of the long pulse, 
and the medium case is t 0 = (D 2 — Di)/4. From Figure 8, we find that when the short pulse is located 
at the edge of the long pulse, the degree of concentration is reduced. Thus, we adopt t 0 = (D 2 — D{)/2. 
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Figure 8. C{y\ t 0 ): the integral of the point spread function (PSF). 
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Selection of r It is desirable for S-BOTDR to suppress the synthetic spectrum's fluctuation as much 
as possible while keeping the peak level high. To meet these requirements, we define the signal-to- 
fluctuation ratio (SFR) as: 



SFR = 



Mpeak 



(45) 



peak 



where /x peak and (J^ eak art the mean and the variance of the synthetic spectrum at the peak, respectively. 
As shown in Appendix C, the SFR is represented explicitly as a function of r, and the optimal r, which 
maximizes the SFR, is given by: 



r opt 



(r^) 2 -2TD 1 + 2(1- e- rDl ) 



1/4 



(46) 



(TD 2 y-2TD 2 + 2(l-e- rD ^) 
Numerical examples of r opt are shown in Table 1. We used the optimal r in simulations and experiments. 

Table 1. Numerical examples of r opt . 



D x [ns] 




D 2 (ns) 




30 


40 


50 


60 


1 


0.0896 


0.0750 


0.0656 


0.0589 


2 


0.1499 


0.1254 


0.1096 


0.0985 


3 


0.2020 


0.1690 


0.1477 


0.1327 


4 


0.2493 


0.2085 


0.1823 


0.1638 


5 


0.2930 


0.2451 


0.2143 


0.1925 


6 


0.3342 


0.2795 


0.2444 


0.2195 


7 


0.3731 


0.3121 


0.2728 


0.2451 


8 


0.4102 


0.3431 


0.3000 


0.2695 


9 


0.4458 


0.3728 


0.3260 


0.2929 


10 


0.4800 


0.4014 


0.3510 


0.3153 



The SFR of S-BOTDR is obtained by substituting (46) into (76) in Appendix C. The SFR values for 
various Di and D 2 are plotted in Figure 9. From this figure, we find that SFR is nearly proportional to 
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Di for fixed D 2 . On the other hand, BFS estimation error variance is inversely proportional to SFR [20]. 
Therefore, the BFS estimation error variance is nearly proportional to 1/D\. 

Figure 9. Signal-to-fluctuation ratio for various D\ and D 2 . 
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4. Evaluation by Simulation 

Numerical simulations were performed to verify the S-BOTDR and compare it with the conventional 
BOTDR. The condition of the optical fiber is shown in Figure 10. In the figure, an assumed BFS is 
shown in (a) and the corresponding Lorentzian spectrum is shown in (b). Simulation results are shown 
in Figures 1 1-13. The parameter values for processing are as follows: 

• Conventional BOTDR (Figures 11 and 12) 
Pulse width is D = 10 and 1 ns. 

• S-BOTDR (Figure 13) 

Pulse widths are D\ = 1 ns for a short pulse and D 2 = 50 ns for a long pulse. The amplitude ratio 
is r = 0.066. The short pulse is located at the center of the long pulse. 

All results were obtained by summing 2 12 trials. For conventional BOTDR with a 10-ns pulse in 
Figure 11, the linewidth of the spectrum was small and smooth; however, spatial resolution was inferior 
and could not detect any interval less than 50 cm. For the conventional BOTDR with a 1-ns pulse 
in Figure 12, because the linewidth of the spectrum was quite large (about 1 GHz), the estimation 
accuracy of the BFS was inferior. In the case of the S-BOTDR in Figure 13, the Brillouin spectrum 
was close to the Lorentzian spectrum; hence, accurate BFS estimates were obtained. Thus, S-BOTDR 
has a high-resolution capability exceeding that of conventional BOTDR. 

The simulation was performed using the exact values of the input parameters. In an actual situation, 
these parameters are, however, subject to errors. As in S-BOTDR four measurements are combined, the 
input light amplitude fluctuation and phase-setting errors become issues of concern. Among them, the 
light amplitude fluctuation becomes negligible by averaging many replicate measurements. The errors in 
the phase differences between short and long elements are, on the other hand, anticipated. We performed 
a series of additional simulations in which the phase differences were subject to random errors. The 
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results clearly demonstrated that they did not affect the spatial resolution, while the BFS estimation error 
increased only slightly, even if phase errors as high as 40 degrees were applied. 

Figure 10. Simulation condition: (a) assumed Brillouin frequency shift (BFS) 
and (b) the Lorentzian spectrum. 



100 MHz 
80 MHz 
60 MHz 
40 MHz 
20 MHz 
0 MHz 
0 



m 1m 2m 3m 4m 5m 



(a) 




v [MHz] 



(b) 



Figure 11. Simulation results of conventional BOTDR with D = 10 ns: (a) the 
Brillouin spectrum and (b) the estimate of BFS. 
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Figure 12. Simulation results of conventional BOTDR with D 
Brillouin spectrum and (b) the estimate of BFS. 
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Figure 13. Simulation results of S-BOTDR with D x = 1 ns and D 2 
(a) the composite Brillouin spectrum and (b) the estimate of BFS. 
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5. Experimental Verification 

We performed BOTDR and S-BOTDR experiments using a test fiber. We constructed a test fiber by 
connecting two kinds of fibers: A and B (Figure 14a). The BFS between fibers A and B differs, and 
the difference corresponds to a strain of 1,200 \ie. In the experiments, fibers A and B are considered 
strained fibers and strain-free fibers, respectively. The replicate measurement numbers were 2 16 in all 
cases. The experimental result for BOTDR is shown in Figure 15. The pulse width of the pump light 
for BOTDR was 5 ns. Strain estimates and the Brillouin spectrum at the 5-m strain interval are plotted 
in Figure 15a,b, respectively. The linewidth of the Brillouin spectrum exceeded 360 MHz (Figure 15b). 
Although intervals larger than or equal to one meter were accurately detected, the strain estimates of the 
other intervals differed from the true values (Figure 15a). 

Figure 14. The fiber under testing used in experiment: (a) the fiber setup and (b) the 
equivalent strain. 
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Figure 15. Experimental BOTDR results: (a) strain estimates and (b) the Brillouin spectrum 
at the 5-m strain interval. 
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We performed S-BOTDR experiments using the following combination of short and long elements 
for the pump light and the low-pass filter: 

Case 1 Di = 7.8 ns, D 2 = 60 ns, r = 0.26 
Case 2 D 1 = 2.7 ns, D 2 = 50 ns, r = 0.14 
Case 3 D 1 = 1.6 ns, D 2 = 50 ns, r = 0.093 

Here, the values of r were determined to be the optimal values. The strain estimates in each case are 
plotted in Figure 16a-c. We attained a spatial resolution of 50 cm in Case 1, 20 cm in Case 2 and 
10 cm in Case 3. Figures 16a-c also indicate that the fluctuation of estimates increases as D\ becomes 
small. This is because the estimation error variance is nearly proportional to 1/Di, as is described in 
Section 3.4. Therefore, some compromise is required to select small D\. The synthetic spectrums at the 
5-m strain interval are shown in Figure 16d,e. In all cases, their linewidths were less than 40 MHz. This 
is because the linewidth of the synthetic spectrums is determined by the reciprocal of the long pulse's 
length. These experimental results prove that we attained 10-cm spatial resolution that exceeds the limit 
of conventional BOTDR. 

Figure 16. Experimental S-BOTDR results: (a-c) Strain estimates of Cases 1 to 3; (d-f) 
composite Brillouin spectrums at the 5-m strain interval of Cases 1 to 3. 
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6. Conclusions 

We showed that conventional BOTDR has, in principle, a resolution limit and proposed a synthetic 
spectrum approach, referred to as S-BOTDR, to overcome this limit. S-BOTDR uses composite pump 
lights and low-pass filters that are composed of short and long elements with different phases. Four 
BOTDR measurements are performed using four specific pairs of phases. By combining four Brillouin 
spectrums with specific weights, an ideal spectrum is obtained. The spectrum is close to the Lorentzian 
spectrum, and the resolution limit, in principle, disappears. We experimentally evaluated S-BOTDR and 
proved that we attained a 10-cm resolution that exceeds the resolution limit of conventional BOTDR. 
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Appendix A. Derivation of (18) 

In this appendix, we derive Equation (18), which is the convolutional representation of the expectation 
of the Brillouin spectrum. Substituting Equation (15) into Equation (16), we have: 




(47) 



whose expectation is calculated using Equation (14), such that: 




(48) 
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We assume that functions that contain ub(z) vanish outside the interval [0, L f ]. The integral interval can 
then be replaced by [—00,00], and we have: 

.2 S ) POO POO POO 



EV(t,u) = V ~^f T T T f(r-s)h(t-r) 

J —OO J —OO J —OO 



OO J —OO J —OO 

x/* (r' - s) h*(t - T^e-^-^^^-^e-^-^dsdrdT' 

/OO POO POO 
/ / f(r-s)h(t-r) 
-OO J —OO J —OO 

xf*(r- t" - s) h*(t -r + T'yT^e- 2n ^-^ s / 2 ^ T ''e- r \ T ''\dsdT'' 

/OO POO 
e -2W / {[f{j »_ s)h{ t- T ")] 
-00 J —00 

* [f\-r" - s)h*(t + t")}} ^e^^^^e-^^l^rf^ (49) 

where we set 7 R = v g ^\Q /2T 2 . Since the Fourier transform of a product is the convolution of the 
corresponding transforms, we have: 

/OO 
{J>[/(r" - s)h(t - t")] FAn-r" - s)h*(t + t")}} * L(s, v)ds 
-OO 

/OO 
\T T „[f(T"-s)h(t-T")}f*L(s,v)ds 
-00 

/OO 
\^ T [f(T)h(t-s-T)}\ 2 *L(s,iy)ds 
-OO 

= lK\^r[f{r)h{t-r)]\ 2 ^ L(t,v) (50) 

where L(t, v) is defined by Equation (19) and is the Fourier transform of (r/2) exp(27riu B (v g t/2)r — 
r|r|) with respect to r. Moreover, defining v) by Equation (20) yields Equation (18). 

Appendix B. Uncertainty Relation of PSF 

The frequency integral of the PSF Equation (20) gives a temporal PSF as follows: 

/OO 
*(t,v)dv=\f(t)\ 2 *\h(t)\ 2 (51) 
-OO 

Center /x T and radius A T of \Pt(*) are obtained as: 

1 f°° 

= jf J t* T {t)dt = Hf + n h (52) 
1 r°° 

A " T = NJ (t ~ ^T(t)dt = a 2 + a 2 h (53) 

where: 

POO 

M = V T (t)dt (54) 
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is a normalizing constant and: 



t\f(t)\ 2 dt 

2 , {t-Hf) 2 \f(t)\ 2 dt 



•oo 
oo 



(55) 
(56) 



Here, || • || denotes a L 2 norm, and ji h and cr^ are defined in the same way. Similarly, a frequency PSF is 
obtained by a time integral as follows: 



(57) 



where / and h are the Fourier images of / and h, respectively. A center and a radius of ^(V) are also 
obtained as: 



A 2 2,2 
' J • n F = a f + a ' h 



(58) 



where fj,;, ny at and o~ h are defined the same as /if and a/. As is well known in Fourier analysis, the 
following uncertainty relations hold: 



1 1 



the use of which, we have: 



A 2 A 2 



aJ + (T2)((7j + 4)> 



47T 



+ 



(59) 



(4tt) 2 (4tt)2 U 2 a 2 



(60) 



Since aj/al + o\ja^ in the last term takes the minimum value of two, when Of = a h , we obtain: 



A 2 A 2 > 



(2tt) = 



(61) 



Appendix C. Signal-to-Fluctuation Ratio of S-BOTDR 



The synthetic spectrum of S-BOTDR Vs(t,v) is the weighted sum of elements V^\t,u),j = 
1,2,3,4, as shown in Equation (32). Since all elements are statistically independent and \cj\ = 1, 
the variance of Vs(t, v) becomes the sum of the variances of V^\t, v), j = 1, 2, 3, 4. Moreover, since 
each V^\t, v) follows an exponential distribution, we have: 



E [V s (t, u) - EV s (t, v)] 2 = J2 ( EvU) (^ v))' 



(62) 



We now arbitrarily fix time t and obtain the fluctuation at the peak of the spectrum. For simplicity, we 
assume that the BFS is constant and u B (v g t/2) = 0. The peak position is then given by v = 0. The 
expectation of all elements of the synthetic spectrum in S-BOTDR at v = 0 is then given by: 



EV (j) (t, 0) = 7r / \L(t, -v) * ¥ j) (t, v 



dv 



(63) 



Sensors 2014, 14 

Using the assumption for the BFS, the Lorentzian spectrum can be written as: 

r 2 



L(t,v) = L 0 (u] 



def 



T 2 + (27H/) 2 

which, substituting into Equation (63), yields: 

/CO POO 
L 0 (u) / ¥ j) (t,v)dtdv 
-oo J — oo 

where we abbreviate the arguments, because EV^ (t, 0) does not depend on t. Noting that: 



L 0 (v) = T r 



-r\ T \ 



and using Plancherel's theorem and Equation (34), we obtain: 



EV {j) 



¥ j \t,v)dt 



dr 



7r 
7R r 



OO p 

i e -r|r| 



[h (j \t)h {j) *{t + T)}dt / [f {j \t)f( j) *(t-T)}dtdr 



/CO 
e" r|T| (K n + re~^K 12 + re i0 'K 21 + r 2 K 22 ) 
-CO 



(K n + re~^K 12 + re^K 21 + r 2 K 22 ) dr 



where: 



K M = 



fk(t)fi(t - r)dt = / h k (t)hi(t + r)dt, fc, I = 1, 2 



Substituting the values of Equation (38) into Equation (67), we have: 

^a) = 2|£( ao + aij . r 2 + a2 /) 

where r is the amplitude ratio of the short and long pulses appearing in Equation (23), and: 



a 0 
a 2 



/CO 
e-^K^dr 
-CO 

pco 

2 / e" r l T l {K n K 22 + cos(% + 4>j)K 2 l2 - K 12 K 21 )dr 



e- T ^K 2 22 dr 



From the relation of Equation (62), cr 2 e3k , the variance of the spectrum at the peak, is represented 



3=1 



into which substituting Equation (63) gives: 



peak 



(f 



4(a 0 + air 2 + a 2 r 4 ) 2 + ^(a 2 ^- - a\)i 

3=1 
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where we set a\ = Yl)=i a i,j/4- 

On the other hand, /i pea k, the mean of the spectrum at the peak, is approximated by Equation (43) as: 

Mpeak = max EV s (t, v) « 8r 2 D 2 ^ R (75) 
Substituting Equations (75) and (74) into Equation (45), we have: 

SFR = (76) 

2 1 4 



Appendix D. Optimal Amplitude Ratio between Short and Long Pulses 

The maximization of the SFR of Equation (76) is equivalent to the minimization of: 

J = ^ + ai + a 2 r 2 (77) 
The value of r that minimizes J is given by: 



'opt , i 
\ a 2/ 



' (78) 
which maximizes SFR. From Equations (22), (23) and (68), 

Kn(r) = ( Cl " |T| ' |T| £ A (79) 
I 0, otherwise 

K 22 (r) = I D2 " |T| ' |r| ^ 2 (80) 
I 0, otherwise 

holds, and by substituting them into Equations (70) and (72), we have: 

a 0 = ^ [(TD 1 ) 2 -2TD 1 + 2(l-e- rD ^)} (81) 
a 2 = ^[(TD 2 ) 2 -2TD 2 + 2(l-e- rD ^)} (82) 

_ (r^-2rA + 2 (i-e-^) \ 1/4 
opt 1 (r^ 2 ) 2 -2r^ 2 + 2(i-e- r ^; 1 • 



Therefore, 



becomes the optimal value of r. 
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